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Abstract 



The properties of strange hadronic matter are studied in the context of the 
modified quark-meson coupling model using two substantially different sets of 
hyperon-hyperon (YY) interactions. The first set is based on the Nijmegen 
hard core potential model D with slightly attractive YY interactions. The 
second potential set is based on the recent SU(3) extension of the Nijmegen 
soft-core potential NSC97 with strongly attractive YY interactions which may 
allow for deeply bound hypernuclear matter. The results show that, for the 
first potential set, the E hyperon does not appear at all in the bulk at any 
baryon density and for all strangeness fractions. The binding energy curves 
of the resulting NA3 system vary smoothly with density and the system is 
stable (or metastable if we include the weak force). However, the situation is 
drastically changed when using the second set where the £ hyperons appear in 
the system at large baryon densities above a critical strangeness fraction. We 
find strange hadronic matter undergoes a first order phase transition from 
a iVAH system to a iVEH for strangeness fractions fs > 1.2 and baryonic 
densities exceeding twice ordinary nuclear matter density. Furthermore, it is 
found that the system built of iVSH is deeply bound. This phase transition 
affects significantly the equation of state which becomes much softer and a 
substantial drop in energy density and pressure are detected as the phase 
transition takes place. 

PACS:21.65.+f, 24.85. +p, 12.39Ba 
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I. INTRODUCTION 



Extremely dense strange quark matter with a strangeness fraction of order one and charge 
neutrality has been suggested to be the absolutely stable form of matter at high densities 
(or at least to be metastable because of the weak interaction) [0-0] . Finite chunks of strange 
quark matter with large strangeness fractions, the so-called strangelets, are predicted to be 
more stable than normal nuclei [§,|5|-0]. Recent theoretical and experimental searches for 
strange quark matter can be found in Refs. 

On the other hand, metastable strange systems with strangeness fractions of order one 
and charge neutrality might also exist in the hadronic phase at moderate values of den- 
sity, between two and three times nuclear matter density. The properties of metastable 
exotic multihypernuclear objects with A and H hyperons reveal quite similar features as the 
strangelets proposed as a unique signature for quark-gluon plasma formation in heavy ion 
collisions [ITO]. The equilibrium between the quark and hadronic phases has been studied in 



ITJ . Strange hadronic matter in bulk has been discussed by Glendenning fT2|j . 

It is found that in an extended mean field theory, a large class of bound multistrange 
objects formed from combinations of (N, A, S) baryons are stable against strong decay. The 



presence of filled A orbitals blocks the strong decay SiV — > AA ||13|| . The maximal binding 
energy per baryon of Eb/A pa — 21MeV occurs at a strangeness fraction or strangeness per 
baryon f s pa 1.0-1.2, charge per baryon f q pa -0.1-0.0 and baryon density 2.5-3 times that of 



ordinary nuclei [13 . It is comparable to that of hypothetically stable strange quark matter 
("strangelets"), which has a binding energy per particle Eb/A pa -10 to -20 MeV. 

The predicted phenomenon of metastability of strange hadronic matter and the actual 
values of the binding energy depend specifically on the partly unknown hyperon potentials 
assumed in dense matter. Some studies used basically Brueckner-Hartree-Fock (BHF) calcu- 
lations with different Nijmegen soft core potentials [|l5| , |T6l . Other studies |T^,|T^,|T8[ extend 
relativistic mean field theory (RMF) from ordinary nuclei (fs = 0) to multistrange nuclei 
with (fs 7^ 0) with the attractive YY interaction of the Nijmegen potential Model D [[HI . 

Recently, Stoks and Lee 0] have studied strange hadronic matter using BHF theory 
and G matrices for coupled baryon channels using an 577(3) extension [pi|-|2"3"|] of the Ni- 
jmegen soft core NSC97 potential from the S =0,-1 sector into the unexplored S =-2,-3,-4 
sector. They have shown that NAS systems are only loosely bound and that charge neu- 
tral strangeness-rich hadronic systems are unlikely to exist in nature. Unfortunately their 
procedure is not self-consistent due to the constraint of equal hyperon fractions that they 
impose and hence does not give the minimum energy ]21| . 



Schaffner-Bielich and Gal [g4j have carried out unconstrained RMF calculations and 
have found larger binding energies. For small strangeness fractions fs < 1, strange hadronic 
matter is mainly composed of NAS and the calculated binding energy closely follows that 
calculated by using the hyperon potentials of the earlier versions of the soft core Nijmegen 
potential. For larger strangeness fractions fs > 1, the calculated binding energy increases 
substantially due to a phase transition into N'EE dominated matter. 

In this work we use the quark-meson coupling Model (QMC) to study strange hadronic 
matter. The QMC model uses the quark degrees of freedom explicitly by coupling the scalar 
a and vector uj mesonic fields directly to the up and down quarks inside the hadrons which 
are treated as non-overlapping MIT bags. It was first proposed by Guichon [p5| and has been 
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used in various nuclear calculations p6|j27[|. More recently it has been modified to include 



a medium-dependent bag parameter coupled to the scalar a field [ 28j29| . This coupling 



is motivated by invoking the nontopological soliton model for the nucleon [BO . Quark 
deconfinement in hot nuclear matter and the phase transition from the hadronic phase to 
the quark-gluon plasma have been studied in this model PJ],P3] • This modified quark-meson 



coupling model (MQMC) f2q]29|| has been extended to the study of strange hadronic matter 



by introducing additional scalar a*(ss) and vector 4>(ss) mesonic fields which couple only 
to the strange quark. This extended model has been used to study the properties of hot 



hypernuclear matter and neutron star matter 

We assume strange hadronic matter to consist of the baryon members of the SU(3) 
octet and decuplet. The octet is comprised of the nucleons and the A, E and 2 hyperons 
while the decuplet contains the A, E, 2 and Q baryons. The decuplet baryons are however 
not found to have a significant contribution in the current calculations at zero temperature 
and are henceforth dropped from our formalism. We studied this system with two sets of 
hyperon-hyperon (YY) potentials, the first set is determined from the Nijmegen hard-core 
potential Model D and the second set corresponds to the potentials obtained in a recent 
SU(3) extension of the Nijmegen soft-core potential Model NSC97. The differences between 
these sets are essentially attributable to the extremely attractive EE and 22 interactions 
in the second set which allow the possibility of deeply bound nuclear matter with hyperons. 
We also consider these two potential models in the context of the QMC model with medium 
independent bag constant. 

The current paper is organized as follows. In section II we describe the details of the 
MQMC model, in section III we show how we fit the various parameters in the model and 
in section IV we present and discuss our results. Finally we summarize our conclusions in 
section V. 



II. THE MQMC MODEL 

In the QMC model, the quark field ip q (f,t) inside a bag of radius Ri representing a 
baryon of species i satisfies the Dirac equation 

V,(r,t) = 0, (1) 

where the quark is assumed to couple directly to the scalar and vector meson fields and m q 
is the current mass of a quark of flavor q. The current quark masses are taken, for the up 
and down flavor quarks, to be m u = = while for the strange flavor m s = 150MeV 
Inclusion of small current quark masses for the nonstrange flavors or other values for the 
strange flavor leads only to small numerical refinements of the present results. In the mean 
field approximation the meson fields are treated classically and the space like components 
of the vector fields vanish for infinite systems due to rotational invariance. As a result 
cj m 7 m =< ujq > 7° = co>7° and =< 0o > 7° = 07° • The nonstrange (up and down) 

flavor quarks are coupled to the scalar er(550) and vector u(780) mesons while the strange 
flavor quarks are coupled to a* (975) and 0(1020). 

For a given value of the bag radius Ri for baryon species i and the scalar fields a and a*, 
the quark momentum x % is determined by the boundary condition of confinement which, for 
quarks of flavor q in a spherical bag, reduces to jo(x q ) = (3 q ji(x q ), where 



ij^dfj, -m q + {g qa a - ^ayy") + (g qir *cr* - g q ^Y) 
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We have defined the effective quark mass inside the bag as 



m q = m q 



and the effective quark energy is given by 



n\(<T,<T*)/Ri = yJ(x q /Ri) 2 + m q 
The bag energy for baryon species i is given by 
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where term is the zero-point energy of the quarks and Bi(a, a*) is the medium-dependent 
bag parameter. In the simple QMC model, the bag parameter is taken as Bo corresponding 
to its value for a free baryon. The medium effects are taken into account in the MQMC 



model @,^] by coupling the bag parameter to the scalar meson fields. In the present work 
we use the following generalized ansatz for the coupling of the bag parameter to the scalar 
fields 



Biip, a*) 



B Q exp 



, , bag . bag *n 



(6) 



with <?k a ^ and as additional parameters. Here it may be worth mentioning that in 



Ref. [|18j the bag constant is coupled to the nonstrange a scalar field but not to a*. 
The spurious center-of-mass energy is subtracted to obtain the effective baryon mass 



M* 



pi — <T rt 2 
bag fern ' 



where 



< v > 



3 

J2 n <i x V R l 



(7) 



(8) 



The bag radius Ri for baryon species % is obtained through the minimization of the baryon 



mass with respect to the bag radius |27 



dM* 
~dR~ 



0. 



(9) 



The zero-point energy parameters Z\ of Eq. (4) are used to fit the actual masses of the free 
baryons Mi = 939, 1116, 1189 and 1315MeV and are found to be Z i= 2.03, 1.814, 1.629 
and 1.505 for the N,A, S and H hyperons respectively, corresponding to a free baryon bag 
parameter Bq = (188.1) 4 MeV 4 and a free nucleon bag radius Rq = 0.6fm. 
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The total energy density of cold infinite strange hadronic matter at finite baryon density 
p B is given by 

* = E (2^)3 / ^ V^ 2 + -k) + + \m\^ + + \m\^ (fO) 

where the summation % runs over the 8 members of the baryon octet which reduces for 
symmetric hypernuclear matter to 4 species with the spin-isospin degeneracy factors 7, = 
4, 2, 6 and 4 for N, A, £ and S, respectively and where is the step function at the Fermi 
momentum k l F . 

The effective Fermi energy of baryon species % is given by 



e*(k F ) = yjky + M* + g^u + g i<p <f> (11) 

which is equal to the chemical potential \X{ 

Hi = Bins + Sifxs, (12) 

where Bi and Si are the baryon and strangeness quantum numbers, respectively. The two 
independent chemical potentials \ib and jis are obtained from the conservation of the total 
baryon density 

Pb~E B ai J d 3 k6(k F - k), (13) 
and the total strangeness density 

Ps = 7^3 E Sat J d 3 k6(k F - k). (14) 
The vector mean fields are determined by 

9u 



i m i 



and 

<t> = E (16) 

i m 4> 

where g iu> and g i( p are the meson-baryon coupling constants defined in Eqs. (18) and (19) 
below. 

The pressure is the negative of the grand thermodynamic potential density and is given 

by 

P = I E (^3 / d 3 k^6(k F -k) + + - l -m^ - \n\^\ (17) 
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III. FITTING PARAMETERS FOR YY POTENTIALS 



We assume that the a and to mesons couple only to the up and down quarks while a* 
and 4> couple to the strange quark. We thus set g r ^ = g ra * = g sa = g su = where r refers 
to the up and down flavors while s denotes the strange flavor. By assuming the SU(6) 
symmetry of the simple quark model we have the relations g sa * = \pig ra and g S( p = \/2g rL0 . 
The coupling of each baryon species with the vector mesons is calculated by counting the 
constituent quarks 

3 

9iu> =J2dqu>=J2 9™, (18) 

q r 

and 

3 

9ict> = ^9q<f> = ^9sct>- (19) 
q s 

With these assumptions the only free parameters left at our disposal are the quark-meson 

coupling constants g ra and g ruJ and the bag coupling constants g^^, gf^ for the 4 baryon 
species and these parameters are adjusted to fit nuclear properties as well as the extrapolated 
properties of hypernuclear matter. The coupling constants of the scalar and vector mesons 
to the nonstrange quarks are taken as g ra = 1 and g ruJ = 2.705 which together with a bag 

coupling constant = 6.81 yield a binding energy of 16 MeV and a compressibility Ky 1 
of 289 MeV at the normal saturation density po — 0.17fm~ 3 of nuclear matter p9|j3T| . 

Table I summarizes the values used in the current work for the basic quark-meson cou- 
pling constants as well as the two sets of the coupling constants gf^ and in the bag 
parameter Bi(a,a*). These sets are chosen to fit nuclear and hypernuclear properties. The 

parameters g ia are taken to fit the hyperon potentials in nuclear matter: 

U ( A N \p ) = -30MeV, 
U^ N \ Po ) = +30MeV } 
ui N \p ) = -18MeV, 

(20) 

where the hyperon potentials are defined by 

UP = (M* - M { ) + (g iui cu + g i4> 4>). (21) 

However, we make two different choices for the constants gf^- In the first set, referred to 

hereafter as MQMC-I, the medium constants gf^ are adjusted so that the potential of a 
single hyperon embedded in a bath of 5 matter becomes 

t/f (po) = Uf(p ) = -AOMeV (22) 

in accordance with the attractive hyperon-hyperon interaction of the Nijmegen potential 
Model D |T3|Jl8|| . Furthermore, we adopt the approximation U~(p ) ~ U^(p ) to fit the 
medium constants. The resulting (po/2) is about -20 MeV. 
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In MQMC-II, we adopt the YY interactions which occur in the most recent SU(3) 
extension of the Nijmegen soft-core potential Model NSC97 [20-22||. The phenomenology 
in this model departs substantially from that in MQMC-I. In particular, the EE and 22 
interactions are predicted to be highly attractive in some channels, leading to bound states. 

We have adjusted the bag constants where i = A, E, H to reproduce qualitatively the 
same binding energy curves of each hyperon species in its own hyperonic matter B\ as those 
produced by the Model NSC97f |20|j2^| . The resulting binding energy curves are shown in 
Fig.|T]. No binding occurs for A hyperons which are already unbound by 8-10 MeV at the 
rather low density ps = 0.05fm -3 . On the other hand, E matter is deeply bound at -33 
MeV per baryon at p^o which is twice as deep as ordinary nuclear matter, and 5 matter has 
an energy of -23 MeV per baryon at p~ . 

Furthermore, we have also considered the original (unmodified) QMC model where the 
bag constant is considered to be medium-independent and simply takes its free space value 
B . In analogy with the MQMC calculations we consider two sets of parameters. The 
QMC-I set reproduces the YY potentials corresponding to the Nijmegen hard core potential 
model D as in MQMC-I. The scalar coupling constants are taken as g qcr = 5.31 and g scr * = 
V2g ga while the vector ones are taken as g guJ =lA71, 1.868, 5.148 and 1.662 and g S( j,=0.0, 
2.289, 0.0 and 3.017 for N,A,E and S, respectively. The QMC-II set reproduces the YY 
potentials corresponding to the Nijmegen soft-core potential NSC97 as in MQMC-II. The 
scalar coupling constants are taken to be the same as in QMC-I while the vector ones are 
taken as ^=1.471, 2.942, 1.103 and 1.721 and ^=0.0, 4.160, 1.560 and 2.434 for N, A, E 
and 5, respectively. 



IV. RESULTS AND DISCUSSIONS 

The resulting energy per particle curves of each baryon species j in its own matter Ej 
are depicted in Fig. [I] as a function of density for models MQMC-I and II. For nucleons, the 
curves are the same in both models since they fit the same nuclear matter properties. In 
MQMC-I, the interactions for A and E are repulsive with E^=+10 MeV and E^=+5 MeV 
at pb = 0.15fm -3 for A and E hyperons respectively. The A hyperon has a shallow local 
minimum at ps = 0.13fm~ 3 while E£ has an absolute minimum of — 7.5MeV per baryon 
at p ^ 0.29fm" 3 . On the other hand, for MQMC-II, the energy Efc reaches +20 MeV 
already at a density pn = O.lfm -3 . This more repulsive potential can be attributed to the 
very weak underlying AA interaction in the extended NSC97 Model. In contrast, the E 
hyperon is deeply bound with a minimum of —33 MeV per baryon at ps = 0.46fm~ 3 . Also 
H hyperon matter has a strong attractive potential with a minimum of -23 MeV per baryon 
at p = 0.39fm~ 3 . This shows that MQMC-II predicts that the EE and SS interactions are 
highly attractive in some channels, leading to bound states. It is also clear that a mixture 
of E and 5 matter must be very deeply bound unless there is an overwhelmingly repulsive 
interaction between the E and H hyperons. 

Figs. and |3| display the effective masses of the baryons (N, A, E, 5) vs ps for several 
strangeness fractions fg using MQMC-I and II, respectively. In MQMC-I, the effective mass 
for each baryon species decreases monotonically with density. It is seen that as fg increases, 
Mjy increases while M£ and M| decrease. The A hyperon effective mass however is only 
weakly sensitive to variations in fg. In MQMC-II, the situation is rather different, as can 
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be seen from Fig. 0. At first the effective masses of the nucleoli and A hyperon decrease 
smoothly with ps for low strangeness fractions. As fg exceeds 1.4, these effective masses 
jump suddenly to higher values as ps reaches a critical value and then start to decrease 
again monotonically for higher values of ps- The jump or discontinuity in the effective mass 
of the nucleon can be as high as 75 MeV for fg = 1.5 at a critical density of 0.57fm -3 . The 
A hyperon has a jump of around 40 MeV for fg = 1.5 at a critical density ps = 0.57fm -3 . 
The effective masses of the E and 5 hyperons have the opposite behavior with M£ and M| 
decreasing when the net strangeness fraction fg of the system increases. Furthermore they 
suddenly jump to lower values as ps exceeds a critical density for fg > 1.4. It is interesting 
to note that in MQMC-II the various effective baryon masses become independent of fg for 
fs > 1-5 and ps > 0.60fm -3 . This behavior indicates that a phase transition takes place in 
MQMC-II at the critical densities where the discontinuities in the effective mass occur. The 
nature of this transition from the A^AS phase to the A^EH phase will become apparent in 
the following figures. 

Figs. f| and || display the density dependence of the bag radius Ri for the N, A, E and 
H baryons for several strangeness fractions fg using models MQMC-I and II, respectively. 
The thick solid line indicates the limiting bag radius value i?exc = (^/^Pb) 1 ^ 3 where the 
excluded volume occupies all space and the assumption of nonoverlapping bags definitely 
breaks down. It is interesting to note here that the bag radius in MQMC tends to increase 
with the increasing baryon density in contrary to the ordinary QMC models where the bag 
radius tends to decrease slowly or saturate at a constant value. This is because in MQMC 
the bag parameter decreases with density which allows the bag to expand. This problem 
can be cured by introducing quark-quark correlations [32| which however are not included in 
the present calculations. Such correlations tend to make the bags shrink at higher densities 
so that the overlap assumption is not violated and -Rexc wm have values much larger than 
those given by the thick line in Fig. |]. Finally it is observed that in MQMC-I the radius of 
the nucleon decreases while i?s and -R- increase with increasing strangeness fractions while 
R\ is almost unaffected by variations in fg. 

In contrast to MQMC-I where the bag radius increases monotonically with baryon density 
Pb we notice from Fig. |^ that in MQMC-II this variation has sudden jumps or discontinuities 
for fg > 1.4. For example, for a strangeness fraction fg = 1.5, drops from 0.775fm to 
0.725fm when the baryon density reaches ps = 0.55fm -3 . The shrinking in the nucleon bag 
size for fg > 1.4 takes place before the breakdown in the nonoverlap assumption even in 
the current calculations where the quark-quark correlations are not included. Furthermore, 
the phase transition tends to occur at lower and lower densities as fg increases and the 
bag radius becomes independent of fg. The same behavior is found for the radius of the 
A hyperons but the situation is a quite different for the E and H hyperons. R-£ (Re) both 
increase with respect to fg and both jump to a higher value at the critical baryon density 
when the phase transition takes place for fg exceeding 1.4. With fg = 1.5, -Rexc is reached 
at pb = 0.50fm -3 for E and H hyperons indicating that, even with the neglect of quark- 
quark correlations, MQMC-II is valid up to ps = 0.50fm~ 3 for such large finite strangeness 
fractions. 

Fig. |] displays the binding energy per baryon for strange hadronic matter vs the baryon 
density ps for various strangeness fractions fg. It is noted in Fig.|6|(a) that the energy curves 
for MQMC-I have only one global minimum for every value of the strangeness fraction fg. 
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At first, the depth of the minimum increases as fs increases until it reaches a strangeness 
fraction of about 1.3 after which the depth of the minimum decreases again for higher values 
of fs and approaches that of purely H matter for f s = 2. It is found that in MQMC-I 
strange hadronic matter is built out always from a mixture of N, A and H baryons. Hence 
each energy curve has only one global minimum. 

Fig. |](b) indicates that the situation for MQMC-II is drastically different. For low values 
of fs each curve still has a single minimum whose depth increases with the increase in the 
strangeness fraction. However, when fs reaches about 1.3, a second minimum appears at a 
rather high density ps = 0.86fm~ 3 . Moreover, this second energy minimum quickly deepens 
and reaches —39 and —56 MeV for fs = 1.4 and 1.5, respectively. The appearance of this 
second minimum can be attributed to the increasingly dominant presence of the E hyperons 
in the system which seems to undergo a first order phase transition from iVAS and NS to 
VEH strange hadronic matter but at a rather higher baryonic density. For values of the 
strangeness fraction fs greater than 1.5 the potential well depth starts to decrease again as 
the S starts to dominate and the energy curves gradually approach that for pure 5 matter 
for f s - 2. 

Fig. [7] displays the depth of the energy minimum and the corresponding baryonic density 
as functions of the strangeness fraction fs- It is evident that the depth of the minimum in 
MQMC-I increases with the strangeness fraction until the curve has a maximum binding 
energy of 19MeV at fs = 1.2. The binding energy magnitude decreases again as fs increases 
to form pure S matter at fs = 2.0 with a binding energy of 10 MeV. As can be seen from 
Fig.|7|(b) the locations of these minima appear at higher densities as fs increases until it 
reaches fs = 1.2 where it hovers around ps = 0.40fm -3 for fs = 1.2— 1.6. Beyond fs = 1.6, 
the location of the energy minimum moves to lower densities. 

The situation for MQMC-II is rather similar to that of MQMC-I for fs < 1.3, but it is 
essentially different for fs > 1.3. When the strangeness fraction reaches fs = 1.3, a second 
minimum in the binding energy appears at pB = 0.86fm -3 . Furthermore, the first minimum 
disappears at fs — 1.5 and a transition to the second minimum (the VE5 phase) occurs 
due to the reaction NS — > EE. The upper and lower branches of the MQMC-II curve 
in Fig.[F](b) thus refer to two different phases with the discontinuity in the density of the 
system indicating it is a first order transition. The binding energy reaches a maximum of 
56 MeV for fs = 1.5 and ps = 0.64fm~ 3 and then decreases to 25.7 MeV as fs increases to 
fs — 2 where the second minimum disappears since the system is now composed solely of H 
hyperons. We also display in Fig. [I] the results for the QMC-I and QMC-II models which 
are seen to be very similar to those found in MQMC-I with a single stable NAE phase. 
In particular, unlike MQMC-II, QMC-II does not lead to a stable iVES system. Finally 
it is worth mentioning that, due to the shrinking of the bag radius with increasing baryon 
density, the nonoverlap assumption survives up to ps > l.lfm -3 . 

Figs. |] and |9| display the fractional density of each baryon species Pi/ps in strange 
hadronic matter versus the total baryon density pB for several values of fs for Models I 
and II, respectively. It is seen that in MQMC-I the E hyperon does not appear at all while 
the A has a substantial contribution for low densities ps < O.lfm -3 . The A contribution, 
however, drops substantially for higher densities to saturate at a constant value around 
Pa/ Pb = 0.20 — 0.25 for ps > 0.2fm -3 and fs = 0.5 — 1.8. It also disappears for fs = 
(nuclear matter) and fs = 2.0 (pure S matter). The contribution of the nucleons Pn/pb 
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decreases monotonically from one to zero as fs increases from zero to two. For fs > 1 
the nucleons comprise a surprisingly small fraction of the total number of baryons for small 
densities but their contribution increases rapidly and saturates at higher densities. The 
variation of the contribution of the S hyperons with density is similar to that of the nucleons 
but with p=/p_B increasing from zero to one as fs increases from zero to two. 

Fig. |9| shows however that the situation is significantly different in MQMC-II, with the 
S hyperons appearing above some critical density and their contribution rapidly saturating 
at higher densities. The fractional density ps/p^ takes its maximum value of about 50% 
& t fs — 1-5. The A hyperon contribution is large at low baryon densities ps < 0.10 but 
practically disappears for higher densities for fs > 1. The fractional densities for the iV's 
and S's in MQMC-II are similar to those in MQMC-I except for fs > 1.2 where there is a 
depletion at high densities caused by the appearance of the £ hyperons in the system. The 
nucleons even disappear completely at high densities for fs > 1.5. 

Figs. [II] and |TT] display pj/e*(0) vs p# e* is the Fermi energy and pj is the chemical 
potential, eqs. (11) and (12). If the ratio pi/e*(k F = 0) < baryons of species i cannot 
exist in the system. The baryon species first appears with zero fermi momentum kp = 
when this ratio reaches one and then takes a finite fermi momentum as the ratio exceeds 
one. ^From Fig. ITU] it is evident that in MQMC-I the £ hyperons cannot be produced at all 
in strange hadronic matter. Furthermore, iV production is increasingly suppressed while H 
production is enhanced as fs increases. 

The situation is quite different in MQMC-II, Fig. [11], where the E hyperons appear 
in the system for strangeness fractions fs > 1.0. They first appear at p# > 0.85 for 
fs = 1.0 and as fs increases they appear at lower baryon densities. However, the situation is 
opposite for A which is produced abundantly for low strangeness fractions fs < 1.0 and then 
disappears when f s exceeds 1.0. The contribution of the nucleons decreases significantly for 
fs > 1.3 and they disappear completely from the system at fs — 1.5 when the baryon density 
exceeds p# = 0.65fm -3 . Furthermore, the nucleons disappear at lower baryon densities as 
the strangness fraction fs increases. It is thus seen that A" production in MQMC-II is 
influenced drastically by the transition from the A^AS phase to the A^SS phase. On the 
other hand, H hyperon production in MQMC-II is not much different from that in MQMC-I. 

Finally, Fig.[TJ| displays the pressure vs the baryon density p# for Models I and II. MQMC- 
I has the normal van der Waal-type loop at small densities reflecting the gas to liquid phase 
transition. In MQMC-II the pressure has an additional sudden and significant drop at high 
densities ps > 0.55fm~ 3 for strangeness fractions fs > 15. This dramatic drop in the 
pressure is due to the transition from the A^AS phase to the NT,E phase in strange hadronic 
matter and probably provides an experimental clue for it. 



V. SUMMARY AND CONCLUSION 

In the present work, we have studied, within the context of the MQMC model, the 
properties of strange hadronic matter composed of the SU(3) octet baryons N,A, E and H 
using two different models for the hyperon-hyperon interactions. It is interesting to note 
here that the medium dependent bag constant for each baryon species reflect more realistic 
self interaction with the scalar mean fields which are corresponding to the nonlinear terms 
of the scalar mean fields in the normal RMF calculations. 
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MQMC-I is designed to mimic the consequences of the Nijmegen hard-core potential 
model D [p~9| , p~3|JT7f1 . It is nevertheless constrained by A and E nuclear phenomenology, and 
by a few AA hypernuclear species reported to date. It accounts more realistically for the 
attractive AA and AfE interactions, but ignores altogether E hyperons. On the other hand, 
MQMC-II is designed to generate qualitatively similar baryon potentials to those obtained 
in the BHF approximation from the SU(3) extensions of the Nijmegen soft-core NSC97 
potential models. The NSC97 model has been tuned up to reproduce certain characteristics 
of A hypernuclei, particularly its version NSC97f. It yields particularly attractive EE, EE 
and EH interactions. 

In MQMC-I, the E hyperons do not appear at all and the A hyperons have a substantial 
contribution at low densities for medium strangeness fractions fs = 0.5 — 1.8. However 
this contribution drops quickly and saturates for Pb/po > 1 a t a constant fraction density 
slightly larger than 20%. This model produces smooth binding energy curves over the entire 
range of fs and ps and forms (meta-)stable strange hadronic matter in the bulk. 

Due to the exceptionally strong attractive interactions among the E and H hyperons in 
MQMC-II, the A/AH- and A^E-dominated strange hadronic matter that exists for fs < 1.3 is 
replaced by A/EE-dominated matter for fs > 1.3, with binding energies per baryon reaching 
as much as 56 MeV for fs = 1.5. The equation of state becomes much softer and seems to 
allow deeply bound hyperonic matter. The system built out of AEE forms a particularly 
(meta-)stable state of strange hadronic matter in the bulk. When the first order phase 
transition takes place above some critical baryon density and strangeness fraction, the E 
hyperons appear and their contribution saturates rapidly to a constant fraction density that 
takes its maximum value of about 50% at fs = 1.5. It is expected that such a phase 



transition may drastically change the global features of neutron stars [36 



Our results for MQMC-II are qualitatively similar to those of Schaffher-Bielich and Gal 
24| where the same YY potentials of MQMC-II (what they call Model N) are employed 



within RMF theory but with only hadronic degrees of freedom. There are however some 



important quantitative differences. Whereas in Ref. ]24|] it is found that the second energy 
minimum, corresponding to A^ES dominated matter, becomes more stable than the first 
minimum for fs >1, our current results indicate that this takes place at higher strangeness 
fractions fs > 1.3. Furthermore, the binding energy well for the A^EE system has a maxi- 
mum depth of 56 MeV in the present work while Ref. pi finds the largest depth to be 78 



MeV. In both calculations the largest binding for the A^ES system takes place at fs ~ 1.5. 
The fact that the phase transition takes place at larger strangeness fractions may make it 
more difficult to observe experimentally, if it occurs at all. 

Finally, quark-quark correlations and excluded volume effects may become important 
at high baryon densities in particular when the bags start to overlap at baryon densities 
Pb = 2 — 3po- We shall consider these effects on the properties of strange hadronic matter 
in our future work. 
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TABLES 



TABLE I. Fitting parameters 



Fit set 






bag 


bag 


bag 


bag 


bag 


bag 


bag 


9qa 


9qu) 


9 Na 


9ho 


9\„* 


9^a 




9~o 


9~a* 


MQMC-I 


1.0 


2.705 


6.81 


4.22 


5.45 


1.63 


7.26 


2.27 


9.12 


MQMC-II 


1.0 


2.705 


6.81 


4.22 


0.0 


1.63 


10.28 


2.27 


10.17 
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FIGURES 




FIG. 1. Binding energy per nucleon (TV) in nuclear matter, compared to the binding energy 
per hyperon (A, S, S) in its own hyperonic matter. The hyperonic parameters where chosen to fit 
MQMC-I and II as discussed in the text. 
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FIG. 2. The effective mass vs baryon density in MQMC-I for the N, A, S, S baryons. 
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FIG. 3. The effective mass vs baryon density in MQMC-II for the N, A, S, S baryons. 
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FIG. 4. The bag radius vs baryon density in MQMC-I for the N, A, S, E baryons. 



17 





0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Pb (f m" 3 ) 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Pb (fm 3 ) 



FIG. 5. The bag radius vs baryon density in MQMC-II for the N, A, S, S baryons. 
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FIG. 6. The binding energy per baryon of strange hadronic matter for various strangeness 
fractions in (a) MQMC-I and (b) MQMC-II. 
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FIG. 7. Comparison of (a) the depth and (b) the location of the energy minimum of strange 
hadronic matter as functions of the strangeness fraction f s in MQMC-I and MQMC-II as well as 
in QMC-I and QMC-II. 
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FIG. 8. Fractional density Pi/pB vs. pb for each baryon species in MQMC-I with different 
strangeness fractions, (a) N, (b) A, (c) £ and (d) S. 
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FIG. 9. Fractional density Pi/pB vs. pb for each baryon species in MQMC-II with different 
strangeness fractions, (a) N, (b) A, (c) £ and (d) S. 
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FIG. 10. The threshold production ratio fj>i/e*(kF = 0) for each baryon species and for different 
strangeness fractions in MQMC-I. (a) N, (b) A, (c) £ and (d) S. 
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FIG. 11. The threshold production ratio fj>i/e*(kF = 0) for each baryon species and for different 
strangeness fractions in MQMC-II. (a) N, (b) A, (c) S and (d) E. 
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FIG. 12. Pressure of strange hadronic matter vs baryon density for different strangeness frac- 
tions in (a) MQMC-I and (b) MQMC-II. 
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